Two-qubit parametric amplifier: large amplification of weak signals 
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Using numerical simulations, we show that two coupled qubits can amplify a weak signal about 
hundredfold. This can be achieved if the two qubits are biased simultaneously by this weak signal 
and a strong pump signal, both of which having frequencies close to the inter-level transitions in the 
system. The weak signal strongly affects the spectrum generated by the strong pumping drive by 
producing and controlling mixed harmonics with amplitudes of the order of the main harmonic of the 
strong drive. We show that the amplification is robust with respect to noise, with an intensity of the 
order of the weak signal. When deviating from the optimal regime (corresponding to strong qubit 
, coupling and a weak-signal frequency equal to the inter-level transition frequency) the proposed 

amplifier becomes less efficient, but it can still considerably enhance a weak signal (by several tens). 
We therefore propose to use coupled qubits as a combined parametric amplifier and frequency shifter. 
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I. INTRODUCTION 



While the desire of building code-cracking quantum computers [l| remains as elusive as ever since its inception, 
its pursuit has been one of the major contributing factors to the enormous progress achieved in quantum mesoscopic 
physics and quantum nanodevices over the last decade and a half. These efforts have already resulted in the devel- 
opment of a new branch of mesoscopic digital and analogue devices @-Q] and even new types of materials known 
^ as quantum metamaterials [s, 9], allowing to control quantum coherent media. In this article we describe how two 
coupled qubits can be used as a parametric amplifier. 



A. Superconducting qubits 

ON 
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' From the point of view of quantum computing, a qubit is any two-state quantum system which satisfies certain 
control and readout requirements and can therefore be used for the execution of quantum algorithms (see, e.g., 
Ch. 2). Many efforts have been devoted to the study of theoretical quantum information processing and, more recently, 
significant progress has been achieved on experimental aspects of this field. For example, one of the most fascinating 
£\j . results obtained in experimental mesoscopic physics over the past decade has been the realization of several types 
t-H ' of superconducting qubits, which demonstrated many of the qualities required for quantum information processing. 
Moreover, superconducting qubits in quantum electronics have a potentially much wider range of applications than 
just quantum information processing (see, e.g., Refs. 0-0, H, ) , including: single-photon generators, producing 

quantum (squeezed or Fock) states of the electromagnetic field, quantum transmission lines, quantum amplifiers, 
etc. Various types of superconducting qubits have been produced, including charge, flux, and phase qubits. These 
- - -' use different physical mechanisms to control their states and store information. For example, in a charge qubit, the 
state |1) has one extra Cooper pair, compared to the state |0), while for a flux qubit, the two logical states differ 
by a fraction of the magnetic flux quantum. Superconducting qubits are mesoscopic (e.g., their working quantum 
states may differ by dozens of millions of single-particle states) and scalable (i.e., it is possible to link these together). 
Moreover, their fabrication, control, and readout require techniques already well developed in solid state electronics. 

Several superconductin g qu bit designs of various degrees of complexity and performance have been realized (reviewed 
in, e.g., Refs. i 1 1 EMI- The so-called persistent-current flux qubit [13L fl4j combines relative simplicity with 
decent decoherence times and scalability. It consists of a small superconducting loop (approximately 10 fim across) 
interrupted by three Josephson junctions. The flux quantization condition is 

4>i + 4>2 + h + 27r$/$o = 27m, (1) 

where 4>j is the phase difference across the jth junction, $ is the total magnetic flux through the loop, $o = h/2e, 
and n is an integer. This allows to eliminate one of the phases (e.g., 03 ). Due to the small self- inductance of the 
loop, the difference between the flux $ and the external flux <$> x , as well as the magnetic energy of the system, can 
be neglected. On the other hand, one must include the contribution to the energy from the charges on the Josephson 
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junctions, Q^/2Cj, and the Josephson energy —EjCOS(/)j, where Cj and Ej is the capacitance and the Josephson 
energy of the jth junction, respectively. Up to numerical factors, the junction charges are the momenta canonically 
conjugate to the Josephson phase differences (see, e.g., [la ]. §2.3). Introducing variables (j)± = ((f>i ± 02)/2, it is 
straightforward to show that the classical Hamilton function of the system (assuming that two junctions are identical, 
C\ = Ci — C, Ex—Ei— E, and that C 3 = aC, E 3 = aE) is given by 
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where Ti± are the corresponding momenta and M± are determined by the junction capacitances. The potential energy 
term in the interval is periodic and, for the proper choice of a and with <f> x « $0/2, it contains a double well, the 
minima of which are almost degenerate and correspond to the current flowing (counter)clockwise around the loop. 
These states are chosen as physical qubit states. Transitions between them are enabled by quantum tunneling through 
the barrier. The quantum-mechanical Hamiltonian of the persistent-current qubit is obtained from ([2|) by using the 
relation N = —id/d<f> between the superconducting phase and the number of extra Cooper pairs (proportional to the 
charge on the Josephson junction). Its detailed analysis is given in [13lll4|. Truncating the Hilbert space of the system 
to the two lowest states, which can be done due to the strong anharmonicity of the potential in its Hamiltonian 
can be reduced to the standard pseudospin form, 
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Here e is proportional to the external flux through the qubit, and A is determined by the tunneling matrix element 
between the potential minima. Here we note that persistent-current flux qubits have decoherence times in excess 
of 10 [is at the operating frequency ~ 1 GHz, allow successful coherent coupling of several qubits, and show steady 
improvements. Therefore, these superconducting qubits are promising new elements for versatile quantum circuits. 



B. Parametric amplifiers 



A parametric amplifier is based on the idea of a parametric resonance occurring for a linear oscillator with parameters 
oscillating in time. Such an amplifier (l6l | is implemented as a mixer, where an input weak signal is mixed with a 
strong local oscillator signal producing the strong output. 

In recent years, several new mesoscopic systems have been proposed and implemented as parametric amplifiers. 
These systems include small molecules in intense laser fields |l7|, p olaritons in semiconductor microcavities [l8j . 
current- voltage oscillations in SQUIDs [l9| and Josephson junctions [20M22I . Another proposal uses active and tuneable 
metamaterials 23] to amplify weak signals. Most of these proposals wish to develop a very compact parametric 
amplifier, which can even demonstrate quantum amplification [24]. This indicates that quantum qubit systems (and, 
in particular, superconducting qubit systems) could be very promising candidates for mesoscopic parametric amplifiers. 



C. Outline of results 



The main results of this work can be summarized as follows. 

• Consider an ac drive Asm(uit) applied to a four-level quantum system (e.g., a couple of flux qubits (e.g..[26l.[27t) 
placed in an ac magnetic field, see Fig. [1} with aj pump = lo in resonance with a transition between a pair of its 
energy levels. Even though the evolution of the density matrix is described by a linear equation, the spectrum 
of the density matrix elements (Fig. [2]) shows a sequence of peaks corresponding to different harmonics of the 
ac drive (which is commonly thought to occur only in nonlinear systems). More interestingly, the strongest 
oscillations are not at the main frequency w, but at some of its harmonics, forming a hierarchy of resonances, 
which is unusual for nonlinear systems. The explanation of this phenomenon lies in the structure of the master 
equations' set for the density matrix elements, in which the external signal enters multiplicatively rather than 
additively. 

• Such a hierarchy of parametric resonances makes this two qubit system (e.g., two coupled flux qubit) an efficient 
parametric amplifier and frequency shifter, especially if the weak signal has its frequency w woa k = w close to 
another inter- level transition of the system. The latter can be achieved by tuning the qubit coupling (e.g., like 
in Ref. H3, see also Fig. [T] below). In this case, a weak signal ecoswt generates a combination of harmonics 
mu + kd>. The amplitudes of these harmonics increase with e, and can reach a height of the order of the main 



3 



harmonic of the drive, even for very small e. Thus, the weak signal can be amplified by a factor of up to several 
hundred. 

• The amplification effects (Fig. |3|) are not suppressed by a realistic amount of decoherence (dephasing and 
relaxation) in the system. Thus, e.g., currently available superconducting flux qubits with relatively short 
coherence times can be used as nanoscale, coherent amplifiers of weak signals in the frequency range of several 
hundred MHz. 

• Noise (Fig. HI of the order of the signal cannot suppress signal amplification (Fig. [5]), which is also robust (but, 
of course, weaker) with respect to both: changing the frequency of the weak signal (Fig. |5^) and the parameters 
of the coupled-qubit system (Fig. [B^,). 



II. QUANTUM AMPLIFICATION WITH QUBITS 

One of the challenging tasks for which superconducting qubits seem to be well suited is the amplification of a 
weak signal, a crucially important tool for both technological and scientific applications. For this goal, different types 
of linear or nonlinear resonance devices are commonly used (e.g., Ref. [l6l |25[ ). The problem of signal amplification 
becomes very difficult at the nano-scale. Here we demonstrate that two coupled qubits can be employed as a parametric 
amplifier 16] based on the effect of the parametric resonance between a weak signal and quantum oscillations between 
the quantum levels of the system, driven by an external ac signal (the pump signal). While the actual realization 
details of the qubits is irrelevant for the mathematics of the problem, we stress that the implementation of the proposed 
scheme is ideally suited for mesoscopic superconducting qubits because these are very controllable and versatile. 

Two coupled qubits can be described by the Hamiltonian: 

H = -\Y,[^ a i + ^ W <} + 9 °l °l (4) 

.7=1,2 

where a J z and a x are Pauli matrices corresponding to either the first (j — 1) or the second (j — 2) qubit; the eigenstates 
of a\ are the basis states of the jth qubit at zero coupling. 

To be more specific, let us consider two coupled nominally-identical superconducting persistent-current flux qubits 
(e.g., Ref. 13), where each of the latter consist of a superconducting loop interrupted by three Josephson junctions, 
along which persistent-currents can circulate controlled by applied magnetic fluxes. When two such loops are placed 
next to each other, so that they feel each other's magnetic fields, this situation naturally produces the "antiferromag- 
netic" coupling represented in Eq. @ by the a x -a x term, with g > (see, e.g., Ref. [26|). More elaborate designs can 
produce a tunable coupling (see, e.g., Refs. [27l [28h . with the amplitude and sign of g controlled externally by the 
magnetic flux, $ COU pi, through the coupler loop [27| (see Fig. [T|). 

The state of each qubit is controlled by the magnetic flux <$>g' = fe*\t) $o through it, where $o = h/2e is the flux 
quantum. In the vicinity of /J = = 1/2, the ground state of each qubit is a symmetric superposition of states \L) 
and \R) with, respectively, clock- and counterclockwise circulating superconducting currents of the same magnitude 
I p . In the basis {\L), ® {\L), \R)}^ the two-qubit system can be described by the four-level Hamiltonian Q 

with 6j = Ip^oSfe , here Sfi (t) contains both the pump and the input signals. The tunnelling amplitudes Aj arc 
usually fixed by the fabrication process, but can be tuned if one of the qubit junctions is replaced by two junctions 
in parallel, in a dc SQUID configuration (see, e.g., [29|, H3|). The interaction constant g, as mentioned above, can 
be made tuneable using a coupler loop (Fig. [!}. Note that the density matrix spectrum is directly related to the 
immediately measurable current/voltage spectrum in the tank, which was exploited in [3lj to detect Rabi oscillations 
in a flux qubit. 

For simplicity, we consider two identical qubits; that is, we assume Ai = A2 = A. A direct diagonalization of the 
Hamiltonian (0} leads to the inter-level transition frequencies 

w (3) = ^T7 + . 9 , u^=2g, W 

which are tuneable by changing g. Therefore, a/ 1 -*, a/ 4 ) can be adjusted to a desirable frequency, that was used 
below. 

Let us drive the qubits simultaneously by a control ac pump signal, with frequency w pump = ui and amplitude A, 
as well as a weak input signal with frequency w woa k = uj and amplitude e <C A, to be amplified: 



ej (t) = As'm(ujt) + esm(ut) + V2D£j(t) 



(6) 
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where we also take into account a noise term (with intensity D) which can be due to fluctuations in the signals or an 
environmental noise. As an example, we consider white noise with zero mean: = and (£j(t)£i(t r ) — ^jiS(t — t'), 

where S refers to either the Dirac delta function or the Kronecker delta. The qubit density matrix p can be written as 

P=\ n o6 ^(8)a b 2 . (7) 

a.b— 0,x,y,z 

This is a straightforward generalization of the standard representation of the single-qubit density matrix expression 
using the Bloch vector; the components U a b thus constitute what can be called the Bloch tensor. Then, the master 
equation, 

+ f p, (8) 

can be written down directly [see Eq. (|12l) in the Appendix], using the standard approximation for the dissipation 
operator T via the dephasing (r^i,]^), and relaxation (r 1 !,^) rates, to characterize the intrinsic noise in the 
system. Also, for simplicity, hereafter we assume that the relaxation rates are the same for both identical qubits, 
i.e., r^i = r^2 = r</> and Ti — T2 — T r , and the temperature is low enough, resulting in Zt2 = Zti = 1, where 
Zxj — tanh(Aj/2fcsTj) is the equilibrium value of the ^-component of the Bloch vector. 

In the limit of zero coupling, g = 0, there exists a solution of Eqs. (|12|l with no entanglement between the qubits. 
This solution can be written as a direct product of two independent density matrices expressed through their Bloch 
vectors: 

P 3 = g(! + x i a * + Y 3 u v + Z i a *)- ( 9 ) 
The components of the Bloch tensor n a f, are all zero with the exception of 

(u ox ,ii ov ,n oz ) = (x 1 ,Y 1 ,z 1 y, (u xo ,n vo ,n zo ) = (x 2 ,Y 2 ,z 2 ) , (10) 

which are just the separate qubits' Bloch vector components. If the interaction is nonzero, g ^ 0, the entanglement 
between these qubits makes the components of the Bloch tensor nonzero, and such an entangled state can persist for 
some time even if the interaction is switched off later on. 
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III. SIMULATION RESULTS 



Formally, the set of equations (fT2|) might look complicated but these are just a set of fifteen coupled ordinary 
differential equations which can be easily integrated numerically [32j . Indeed, measuring time in units of 1/A, we 
numerically solved Eqs. (1T21) by the Euler method for two coupled qubits, driven by the field ([6]) with w pump = u = uj^ 
and w W eak = u> = u/ 3 ) [see Eq. ((SJ)]. Our numerical integration produces a set of data for H a ,b(t) which can be analyzed 
using a Fourier transform. Thus, we can numerically evaluate the spectrum of Z\ — n oz , defined as 

S z (u) = {\Z x {u)\), (11) 

and the similarly defined spectrum of X\ — Tio x - If one considers a flux qubit as a particular realization of our 
two-qubit parametric amplifier, the n Z-\ (t) and X\{t) (and, thus, their spectra Sz and Sx) can be directly obtained 
by using impedance measurements [331 ] of the circulating current in the first qubit of the proposed device. These 
spectra are shown in Fig. [2] for g/A = 1. Below we will analyze the influence of a weak signal on Sz and Sx- 

When the weak signal is switched off [Fig. [Ha)], the spectrum of Z\ exhibits several peaks corresponding to the 
harmonics, u m — muj, of the applied drive (A/ A = 15). In contrast to the standard nonlinear response, where the 
spectrum starts from u>\ and has harmonic peaks u) m , whose heights decrease with m, in our two-qubit system, S z 
starts with oj 2 and contains only even harmonics. The peak heights show a surprising non-monotonic dependence on 
m, with the two highest peaks occurring at m = 20 and m = 22. A somewhat similar non-monotonic peak dependence 
on the harmonic number is seen for Sx, the spectrum of the off-diagonal matrix element X\ = H ox - Even though 
this spectrum starts with the main harmonic uji, it contains only odd peaks which heights still show a non-monotonic 
dependence on m [Fig. [2jb)]. The super-indices in ,u^ x \ etc., refer to the level-splitting (intrinsic properties) of 
the qubits, while the lower indices in u)x,u 2 , . . . label the harmonics of the response to the input ac monochromatic 
drive. We measure all frequencies in units of the energy-splitting of independent qubits a/ 2 ' (g = 0), since this provides 
a characteristic frequency in the system which is fixed even if one tunes the inter-level frequencies 
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FIG. 1: Schematic diagram of two flux qubits (persistent current qubits) coupled via a coupler loop. Josephson junctions are 
represented by crosses, and their thickness indicate their Josephson critical currents (drawn not to scale). The qubit states and 
the amplitude and sign of the coupling constant g are controlled by the corresponding magnetic fluxes, /j 1 ' 2 ' and / coup i or (in 
units of $o)- Both the pumping drive and the weak signal Asinut + esinujt are generated by a source in the bottom circuit. 
The left circuit is needed to pick up an output signal Zi(t). The top circuit controls the coupling g. Noise, shown by y/2D^i(t), 
is coupled to each qubit. 



All spectral peak heights are proportional to the drive amplitude A. Thus, a weak signal should produce a spectrum 
with strongly-suppressed peaks. This is consistent with our simulations for zero-drive amplitude A = and weak- 
signal amplitude e/A = 0.1. The resulting spectrum has only one peak, with an amplitude about 150 times lower 
than the highest peak of the Sz spectrum for A/ A = 15, e = [Fig. HJa), the weak peak is shown there as "weak" 
and it is almost invisible]. 

When mixing the strong drive with a weak signal, the resulting spectrum Sz at A/ A = 15, e/A = 0.1 (i.e., 
e/A = 1/150) [Fig.[3Ja)] strongly differs from the simple superposition of the two spectra described above, Sz(A/A = 
15, e = 0) and Sz{A = O.e/A = 0.1). It is remarkable that the weak signal produces a considerable change in the 
spectrum: it generates peaks at combination frequencies, kuj + luj, with integer k and I; their heights are determined by 
the weak signal but can be of the order of the highest peak in Sz(A/A — 15, e = 0). Interestingly, the enhancement 
of combination-frequency harmonics occurs by borrowing some energy from the pumping drive, which own harmonics 
(harmonics existing at e = 0) decay when the weak signal is applied. Indeed, if one compares the spectrum at e = 
shown by red dotted lines, with the spectrum at e/A = 0.1 shown by black solid lines in Fig. [3^, one can see that the 
heights of the harmonics at e = are higher that the corresponding peaks at e/A = 0.1, even though the total energy 
pumped in the system is larger when both the pumping drive and the weak signal are applied. 

The heights of the combination-frequency peaks increase with the weak signal amplitude, e, followed by saturation 
at large values of e. Of course, the heights of the combination-frequency peaks tend to zero when e — > 0. It is clearly 
seen [Fig. [3Jb)] that their height can be approximated by a linear function, /3e, for the weak signal amplitude in the 
range < e/A < 0.005, with an amplification coefficient ft of about 100. This level of amplification is remarkable. 

It is also useful to stress that many peaks associated with different combinations of two frequencies, fcw P ump+'w W cak, 
appear in the spectrum of Sz for different integers k and I. In other words, the spectrum has many combination- 
frequency harmonics with different intensities. This allows to pick up the signal on a frequency which better fits an 
available experimental setup. Thus, the proposed two-qubit amplifier is also a frequency shifter, allowing to shift the 
frequency of a weak signal to a desirable frequency range. Note that it is usually hard to predict which peak should 
be the highest one [see Fig. [3] with highest peak marked as /(e)]. However, by measuring the output signal I± of the 
first mixing harmonics, w wca k ± w pump , usually allows to pick up a strongly amplified signal. For instance, the ratio 
of the highest peak I e to I± is about 1.7 in our simulations shown in Fig. [3^,. Therefore, choosing the harmonics 
w woa k ± Wpump would be a good guide to observe the predicted signal amplification. 
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FIG. 2: (Color online) The spectrum Sz of the Zi matrix element (responsible for the occupation of the excited level of the 
first qubit) is shown for dephasing and relaxation given by IVA = T r /A = 10" 3 (we use the same Td, and F r for all results 
reported in this article), (a) When the weak signal amplitude is equal to zero (e = 0) and when the reduced-drive amplitude 
A/ A = 15, the ac monochromatic drive fed into the two-qubit amplifier is converted into a set of even harmonics ui m = mcj, 
with m = 2,4, 6.... Note the absence of response at u>i — cj (no peak there). For comparison, Sz is shown in blue and indicated 
by the blue arrow for zero-drive amplitude A — and weak signal e/A = 0.1. Only one very small peak is hardly seen, indicated 
by the blue arrow, corresponding to 2ix) wca k. (b) Spectrum Sx of the off-diagonal matrix element X\ (zero signal case). The 
same ac monochromatic drive fed into the amplifier is converted into a set of odd harmonics for Sx'- m = 1, 3, 5, ... . Note that 
(a) and (b) show very unusual non-monotonic spectra. The observed non-monotonicity in the system response can be seen as 
a fingerprint of the qubits, characterizing their dynamical nonlinear response. 
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FIG. 3: (Color online) (a) The amplification of a weak signal e/A = 0.1 by a strong drive A/ A = 15 (i.e., e/A = 1/150) can 
be seen in the spectrum of Sz shown by the black solid lines. The spectrum obtained is not a simple superposition of the two 
spectra Sz{A/ A = 15, e = 0) (also shown here by red vertical dotted peaks) and Sz(A — 0, e/A = 0.1) [see also red and blue 
peaks in Fig. HJa)]. Instead, a combination of harmonics appears for kuj + lui, with integer k and I. The height of these peaks is 
almost proportional to fie, for < e/A < 0.005, and it is strongly enhanced by a factor f3 ~ 100. The enhanced peak is marked 
by an arrow with symbol 1(e). Note that this peak is absent in Fig. [2ja). (b) Normalized output amplitude 1(e) /1(A) (ratio 
of the heights of the highest mixed peak 7(e) to the highest peak 1(A) of the spectrum at e = 0) as a function of the reduced 
weak signal amplitude e/A. This dependence shows an almost-linear increase of the peak height with e for small e/A, followed 
by saturation and even decay at relatively large e/A. 
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FIG. 4: (Color online) Trajectories Zi(t) for different noise levels: VD/e = for (a) and s/15/e = 0.066 for (b). All other 
parameters are the same as in Fig. [3^. As seen from (b), the applied white noise considerably affects the time dependence of 
Z\(t), making trajectories quite noisy with respect to the noiseless situation shown in (a). 
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(a) noise ~ 6.6% 




FIG. 5: (Color online) Spectrum S z (uj) for the same parameters used in Fig. [3] and different noise levels y/D/e = 0.066 for (a) 
and 0.2 for (b). One can clearly see the mixed (pump-signal) peaks for several combination-frequencies i;w pump + Zw W eak, even 
at high noise levels, when the trajectories Z\ (t) shown in Fig. [4] are very noisy. 
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FIG. 6: (Color online) Amplification of a weak signal away from the optimal regime, (a) Spectrum S z for off-resonance 
frequency o; WGa k = 1.113u/ 3 ' (g = 1) of a weak signal shown by red thin solid line. All other parameters of the simulations 
are the same as in Fig. [3] The thick black dotted line corresponds to the simulations with zero signal, thus, highlighting 
the combination-frequency amplified harmonics as red lines with no black dots on top. One can see that the amplification of 
off-resonance signal is weaker, but still remarkable with amplification factor ft about 30. (b) Spectrum of S z for weakly coupled 
qubits g/A = 0.1 and weak signal amplitude e/A = 24, as shown by the red thin solid line; all other parameters are as in Fig. [3] 
and Aj A = 12. The same simulations but with zero signal is shown by a thick black dotted line to highlight the amplified 
mixed peaks. The estimated amplification factor /3 is about 10. 



11 



IV. NOISE INFLUENCE 

Since a weak signal can be considerably amplified by a strong drive in our parametric amplifier, one can wonder if 
an uncontrollable noise could also be amplified making the weak signal indistinguishable. To check this we performed 
simulations at different noise levels chosing all other parameters as in Fig. [31 i.e., A/ A = 15, and e/A = 1/150. A 
noise with intensity of the order of 6.6% with respect to the weak signal (i.e., V~D/e ~ 0.066) has already considerably 
affected the time dependence of the measured signal (see, Fig. [4}, but it does not strongly influence the spectrum of 
Sz(ijj), as seen in Fig. [5k. The reason for this is that the noise contains all frequencies (or at least a broad frequency 
spectrum), thus, its energy pumped in the signal harmonic is relatively small. 

Even stronger noise, \D/e = 0.2, is still not enough to suppress the peaks attributed to the weak signal. Moreover, 
a sort of stochastic resonance (increase of the peak heights with noise) [34| is also seen on Fig. [S}d. Therefore, the 
proposed two-qubit parametric amplifier is robust with respect to noise and, moreover, the noise can even be used to 
further amplify the signal. 

V. AWAY FROM THE OPTIMAL REGIME: ROBUST AMPLIFICATION OF TWO-QUBIT 

AMPLIFIER 

In order to realize the strongest amplification of a weak signal (optimal working regime) , a tuneable coupling must 
be used. In other words, the coupling should be tuned to adjust a level splitting frequency u)^' to the frequency w wea k 
of the applied weak signal. This adjustment cannot be always realized. Thus, the amplification of a weak signal of 
an arbitrary frequency (different from the level splitting) is worth studying. As we expected, the amplification of the 
weak signal on the frequency w W cak = 1.113a/ 3 ) =/= uW is weaker (Fig. [5k.) by a factor of about 3.5 (ratio of the highest 
combination- frequency peaks in Figs. [3] and [6]) with respect to the optimal amplification w wca k = u/% i.e., in the case 
when the signal frequency is equal to the level splitting. Nevertheless, this amplification is still strong enough {(3 is 
about 20-30) allowing to use the proposed amplifier for a weak signal with arbitrary frequency. 

A similar situation occurs if the qubit coupling is not strong enough to reach the optimal regime (see Fig. [S}d) . 
Indeed, the amplification of a weak signal (e/A = 0.5) by the strong drive (A/A = 12) is weaker, but still essential 
(j3 about 10) for the coupling g = 0.1. 

VI. CONCLUSIONS 

The spectrum of two coupled qubits driven by an ac signal with the frequency in resonance with inter-level transitions 
has an unusual structure, with a hierarchy of harmonic peaks with heights non-monotonically dependent on the 
harmonic's number. This peak-height hierarchy is a fingerprint of any two-qubit system and can be used to characterize 
both individual qubit parameters as well as the interqubit coupling. 

Exploiting the analogy between a parametric amplifier and a system of two coupled qubits, we propose a method of 
amplification of a weak signal via its mixing with a strong pump signal applied to the two-qubit system. If both signals 
are relatively close to the inter-level transitions in the four-level quantum system (which can be achieved by tuning 
the qubit coupling), then the amplification coefficient can be of the order of 100. When the weak signal frequency 
is different from the inter-level splitting then the amplification is still strong enough, allowing the proposed amplifier 
to work efficiently both in inter-level resonance and off inter-level resonance regimes. Weakening the qubit coupling 
also suppresses a weaker signal enhancement, thus, requiring strongly-coupled qubits for this remarkable parametric 
amplification. 

We also show that noise, which is of the order of a weak signal, can strongly affect the time dependence of the output 
signal Zi(t), but it modifies much weakly the spectrum Sz- Therefore, the proposed amplifier can work efficiently in 
noisy conditions. This large amplification offers a different way of using multiqubit circuits as parametric amplifiers. 
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VIII. APPENDIX: MASTER EQUATIONS FOR THE DENSITY MATRIX COMPONENTS 

The master equation ([7]) can be explicitly written as follows [34| : 

Iloa; = A 2 Ho y — T^H-Ox 

tl y = -A 2 IL 0x + £ 2 (t)Tl 0z - 2gIL xz - T^Oy 

n 0z = -e 2 {t)H 0v +2gIL xv -T 2 (Ilo z -Z T2 ) 

Ilyo = — AiII^o + ei(t)II z o — 2gU zx — TtpiUyo 
tlzo = ~ei{t)U y o + 2gIl yx -T 1 (n z o- Z T i) 

^xx — A 2 H xy + AiH yx — (r^! + T^^lxx (12) 

H xy = -2gU 0z - A 2 U XX + AiU yy + e 2 (t)U xz - (T^i +T ct , 2 )Il X y 

tl yx = -2gIl z0 -A 1 IL xx + A 2 ILyy + ei(t)IL xz -(T ! j > -i+T ( f >2 )ILy X 

tl xz = 2gTl Q y - e 2 (t)Tl xy + AiU yz - (I^i + T 2 )T1 XZ 

tl zx = 2gn y0 -e 1 (t)U yx + A 2 U Z y-{T^ 2 + r 1 )Il zx 

Jlyy = -A^y ~ A 2 Tl yX + e 2 (t)Tly Z + ei(t)Tl Z y ~ (T^i +T,p 2 )llyy 

ii yz = -A x n xz -e 2 (i)n ra + £l (t)n zz - (r 01 + r 2 )n yz 
ii zy = -A 2 n zx - e^Hyy + e 2 (t)n zz - (Ti + r 02 )n^ 
ti zz = -ei(t)n„, - e 2 (t)n zy - {T t + r 2 ){n zz - Z T1 Z T2 ) 

where the symbols are explained in Eqs. p] [TU)) . In the absence of qubit-qubit coupling, g — 0, the first three 
equations in (JT2J) describe the evolution of the Bloch vector components (Ho x , n 0y , n 0z ) of qubit 1, and the second 
three equations in (fT2"|) describe those (n^o, n y o, II z o) of qubit 2. 
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